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Abstract 

tSJTT 77 on *r *“““ "** ** **" w»*» (Mm mwmojm 

L, ir^ln 7 T7 deacribed in “ s TOrk is responsibie tor combMns hi » h p**>™ 

- and GH. to prodnce a hydrogen Sow that meets certain thermodynamic properties. The 
sire properties are maintained by precise control of the LH2 and GH2 flows. The mixer is 
modeled as a genera multi-flow lumped volume for single constituent fluids using density Id 
ernal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK 

— " 7 s a 7 Ie l00t " U,: ’ featale ° f tbe flUid thermod ^™ Properties. A small- 
gnal (linear) model is developed based on the nonlinear model and simulated as well Pulse 

disturbances are introduced to the valve positions and the quality of the linear model is ascertained 

by comparing its behavior against the nonlinear model simulations. Valve control strategies that 

mu ate an operator-, n-the-Ioop scenario are then explored demonstrating the need for automatic 

are d “d h J l’ ° PtimSl and “““'Output Proportional/Integral controllers 

e designed based on the linear model and applied to the nonlinear model with excellent results 

track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer 

pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve 


1 Introduction 

. The Test and Engineering Directorate at NASA John C. Stennis Space Center (SSC) continues 
ts Sorts to assemble a software ampliation package that captures the static and dynamic charac- 

- " 1CS 01 m ° e ‘ n aS mLUre ibermodynamic systems. The package is foreseen to fulfill the need 
or an accurate and verifiable thermodynamic system simulation environment [3. 4. 5. 6], 

In this work, we build mainly upon a previous report [6] and focus on the point-to-point tracking 
control problem for the LH 2 and GH 2 mixer subsystem and associated flow controllers. The relevant 
components of the mixer subsjrstem are shown in the flow schematic of Figure 1. The LH 2 valve 
com, o.s t e ow of high pressure liquid S 2 from a pressurised tank; the GB 2 valve controls the flow 
o gaseous ife from a set of high pressure bottles; and the output valve controls the flow of H 2 into the 
test article. The primary objective of the mixer and control valves is to regulate the thermodvnamic 
properties of the out flow of H 2 for optimal performance in spite of the steady depletion of the source 
bottles, measurement errors, modelling imperfections, and other uncertainties and perturbations. 


2 Background 

Some application examples of the use of a related software package called EASY/ROCETS 
are presented m [3] including FORTRAN listings of several modules and a brief user’s manual. 

• The report [4] performs a detailed derivation of the model of the RUN-TANK module as 
well as its implementation within EASY/ROCETS. The low pressure and high pressure LOX 
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Figure 1: Diagram of HPH2 Mixer and Control Valves. 


and hydrogen subsystems are simulated with EASY/ROCETS and compared against previous 
computer simulations with good results. In addition, various control strategies were simulated 
for the hign pressure hydrogen subsystem. The objective being the temperature regulation of 

f.b P nntnnf cfroom A 4- L. . t » . . 

lunvever, control aesigns are done in an ad-hoc manner due to 
the lack of a suitable control-oriented mathematical model of the mixer. 


In [5], several upgrades and enhancements to EASY/ROCETS were reported including the 
NIST-12 thermodynamic database. 


Most recently, Victor Marrero, one of the SSC 2001 Summer Undergraduate students, reported 
his work on validating the EASY/ROCETS models of several components of the high-pressure 
hydrogen system with real test data. The results are very encouraging and prompts for 
continued analytical studies such as the one presented in this report. 

• In [6], we developed an analytical model of the mixer and control valves suitable for dy nami c 
analysis and control design in the vicinity of a thermodynamic equilibrium point. An optimal 
linear-quadratic regulator was designed and simulated using Matlab to illustrate the effective- 
ness of such a controller in compensating against perturbations that result in small deviations 
around the desired equilibrium. A simulation was included that indicates how one could study 
which valve or combination of valves is more effective in achieving a desired performance. In 
essence, this is a study on control authority or relegation. 


Results 

Based on. the work by H. Richter [7], the nonlinear mixer equations were simulated in the 
SIMULINK environment. 




• As discussed in [6], there was a need to describe the gas valve using a compressible flow model. 
This is now incorporated in this report and the corresponding changes in the small-signal model 
were derived. The details are in Appendix A. 

. Also based on the work by H. Richter [7], the NIST-12 property data is now incorporated in 
the model. The program MIPROPS was used to generate the necessary tables which are then 
continuously accessed via Matlab functions. 

• My student, Ms. Jamie Granger, as part of her Summer LAMP program with Tulane Uni- 
versity, verified the accuracy of the data in the table and also combined two routines into 
one. 


• The small-signal model developed in -[6] requires the numerical computation of several partial 
derivative terms which in turn depend on the equilibrium point. Since there are no analytical 
expressions that would allow the direct calculation of some of these partial derivatives, a 
Matlab routine was written to compute these terms. 

• Open-loop control strategies were investigated mainly to test that the model does provide 

reasonable results. In the future, we hope to validate the model when real test data becomes 
available. 


. Closed-loop control strategies were investigated including: 1) an optimal multi-input/single- 

output (MISO) PI controller to track the exit flow; and 2) an optimal multi-input/multi-output 

(MIMO) PI controller to f'.ra/'lf j - j__ . i,i ... 

— uxuijci temperature, ana me exit temperature 

simultaneously and in the presence of step perturbations in the valve positions. 


3 Updated Mixer and Valve Models 

The mixer has been modelled as a general multi-flow lumped volume for single constituent fluids 

using density and internal energy as states [3, 6]. The model was independently derived in [7], In 

this report, we refer to the mixers internal energy and density as the “state”, not to be confused 

with the definition of state of a system as described by its thermodynamic properties. Letting 

21 221 ® denote the mker ’ s s P ecific internal energy and density, respectively, the two dynamic 
equations are given by 2 

Zl = yTJ) ~ WnetZi + (jj ^ 


where E net is the net (input minus output) energy into the mixer, W net is the net (input minus 
output) flow, and Q is the total heat transfer rate. The net flow is simply 

Whet — W g + Wi — W e 

2 In [6] the mixer states were its internal energy and its mass (Vz 2 ). 



where the subscripts stand for “gas”, “liquid”, and “exit”, respectively. In this study, only positive 
h ows are co nside red, th at is , t h e ga s and liquid flo ws are into the mixe r a nd the output flow is out 
of the mixer. In essence, normal operation requires that each valve always experience a positive 
pressure difference across it thus disallowing flow reversal. The net input energy term is then riven 


-Bnet 


hgWg -r h[Wi — h v w e 


where h„ h,, h, denote the enthalpy of the GH 2 stream, LH t stream, and volnme block (miner), 
respectively. 

The valve block calculates the incompressible fluid flow value w f given the inlet P { and outlet 
P 0 pressures, the density of the inlet stream p { , and the valve flow coefficient C vf , as follows: 

W f ~ \f Pi Pin (Pi ~ Po) C vf — CZfCyf /Q\ 


v nere p w is the density of water (62.4 Ibm/ft 3 ), and we have assumed that P £ > P Q . As pointed 
out m [6] the model of the gas valve had to be modified to simulate compressible flows. According 
to [7], the following expression models the compressible mass flow through a valve: 


w 


, = ( 1857 x 10~ 2 ^i ] /p/-P i * Cvfi Pq < p. < 2P 0 - 

2.423 x 10 ~ 2 ^T iPi C vf , Pi>2P 0 . 


(4) 


when the mess flow is in Ibmfs. nrp«nrA m nc-j fko a 

l ~7 C ~ ***■ p't-j O A UiXV^ \J. 


density m and the temperature 


T in ° R ■ By desi S n and under normal operation, the properties of the inlet hydrogen gas and 
liquid are such that P t , P g > P 0 = p v an d also P v > P e ; thus, the possibility of reverse flow is not 
considered in this study but can certainly be added if\so -required. 


3.1 Small-Signal Model 

Substituting the expressions for the valve flows and the net internal energy and net flow into 
equations (l)-(2), and treating the valve coefficients as inputs to be regulated, the dynamic model 
of the mixer is a 2-state, 3-input system of nonlinear differential equations of the form 

P ~ Pl(^1; z 2-, C V g , C v l, C ve ) ^ 

^2 ~ p2 ( z l: z 2: C V g, C v l , C ve ) (Qj 

where Fi(-) and p 2 (-) are nonlinear functions of the state and valve coefficients. 

For the remainder of the report, we will denote constant or equilibrium values of any variable by 
an upper bar (•). Given constant values of valve flow coefficients C v = [C vg C vi C ve ] T (superscript 
T denotes transposition) and constant fluid properties, the state of the model 

Internal Energy 
Density 



L ^2 


reaches a constant equilibrium point z. Next, consider perturbing such an equilibrium by small 
signals x(t) and u(t) so that 

z(t) = z + x(t) and C v {t) = C v + u(t) 

where u(t) denotes a small, valve-coefficient correction/regulation signal. Then, a standard lin- 
earization of equations (5)-(6) results in the small signal model 

x = Ax + Bu ^ 

where x (t) is the small perturbation state vector, u(t) is the small perturbation control signal, and 
the two-by-two matrix A and two-by-three matrix B are given by 

dF dF, dF, 

OOyg 9C v l SCje 
dF? dF ? dF-> 
dC-ug BCyl dCve 

where the partial derivatives are evaluated at the equilibrium state z and constant valve flow 
coefficient vector C v . 

Typically, in the control literature an output equation of the form 



' dF_ 

dF ' 


A = 

OF 

dF't 

si; 

dF 

: B = 


dz l 

dz? . 



y = Cx -f Du 


where matrices C and D are appropriately dimensioned is appended to the model (7) to account for 
the measurement of certain variables such as temperature, pressure, or flow. An output of interest 
is the exit flow w e . Using equation (3) and linearizing around a chosen equilibrium point, we obt ain 
the linear approximation 

r da e da e ] 

X' T [ttfijgg. C'vfloW'X A DflrrujUy 


W elin 


dz\ 




eg 


A second output of interest is the mixer pressure P v . Taking this pressure to be a function of the 
mixer internal states, that is, 


Pv 




P olin — 


then, the linearization of this function around the chosen equilibrium gives 

'dP v dP v ] 
oz x dz 2 J eq v 

A third output of interest is the exit temperature T e . Taking this temperature to be a function of 
the mixer internal states, that is. 


T e = fT(zi(i),z 2 (t}) 


then, the linearization of this function around the chosen equilibrium gives 


T e lin — 


'dTe 

_dzi 


dT e 

dz2. 


eq 


X = C te X 


The indicated partial derivative terms are calculated numerically using Matlab routines. 



It is not reasonable to assume that the mixer states be measurable. However, the state postulate 3 
— implies th at the state z(t) =z + x is directly available for control design if at least two independent 
thermodynamic properties were measured. Therefore, 

• assume that real-time measurements of temperature T e and pressure P e are available at the 
exit point. Then, using the thermodynamic table, one finds the exit flow enthalpy h e . 

• Assuming a lossless valve, then enthalpy across the valve is preserved so that h v = h e , where 
K is the mixer (volume) enthalpy. 

Using h v and T v (measured) and the table, one finds the mixer pressure P v . 

• From h v and P v and the table, one finds the mixer states (energy, density). Or, using the 
measured exit flow w e , the known exit pressure P e , the computed P v , and the exit valve 
coefficient, one solves Equation (3) for the density p v . 

The above steps were verified using the program MIP-ROPS. In summary then, it is reasonable to 
■assume that 


= X = Z " f ■ (8) 

is available or that it can be deduced from other real-time measurements. This may be called a 
“static observer' 1 . An alternative is to design a so called dynamic observer to provide a real-time, 
optimal estimate of the states. This is beyond the scope of this report. 

The calculation of the equilibrium points and the elements of thp matri™ A o 
m [6j. me relevant results are compiled in Appendix A of this report. The only difference here is 
that the model of the flow in the gas valve is for compressible fluids. 

The evaluation of is accomplished numerically via a Matlab routine. The evaluation of % 
on the other hand, is done by differentiating the expression 

p 

^ E ~r C units ■ 

P 

where C units = 1 in the SI units or C units = 0.1852 when h, E are in BTU/lbm , P in psi, and p 
in Ibm/ ft 3 . This results in 


dh v 

dzi 

dh v 

dz-? 


= 1 + 


Cunits dP v 


ZO 


= a 


units 


dz\ 
1 dP v 


zo dzo 


Pv 

4 


The relations between pressure P and internal energy or density are essentially linear in the vicinity 
of an equilibrium point. This verifies the work reported in [6], Therefore, now we have the capability 
to compute the small-signal model at any desired equilibrium point. 


3 For a simple compressible substance, any thermodynamic property is at most a function of two other independen 
thermodynamic properties. 



3.2 Purpose of the Small-Signal Model 


e m° e (,) escribes the dynamics of the mixer model in the memity of the eqmSbrmm iinder 
consi era ion. eally, both x(t) and u(t) are zero; in practice, the state z(t) deviates from the 
re z us necessitating a corrective action by u(t). The model is suitable for a variety of 
ana ysis s udies such as stability, controllability, and regnlation/tracking of the equilibrium. The 
stability study is relevant since it reveals an important characteristic of the equilibrium a under no 
con to, that ri, with „(„ = 0. Appropriate action must be taken by the control „(t) in the event 

wheth ““! t °“ e) 01 Wen “ P °° riy teWd ”- The Mti °" °f controllability reveals 

suitable control u(t) exists such that the perturbation state x(t) may be steered to zero 

oreover ,t is possible to study which valve or combination of valves is most effective [6|. Finally 
e equilibrium regnlation/tracking study refers to the actual control design to regulate aft) to zero 
or oo track a desired equilibrium trajectory as a function of time. 

. S ' rateSy tiW ‘ * “ Sed to contrd s ^ tems "*«* “near model changes from one equilibrium 
pomt to the next is known as gain scheduling. Basically, one designs a controller for each small-signal 
mooel m the vanity of each desired equilibrium. An outer-loop controller or scheduler performs 
the decision as to when to switch from one controller to the next. 


3.3 A Numerical Example 

It is to be noted that the thermodynamic tables used by the Matlab routines were generated with 
tne program MIPROFS and therefore contain an inherent error because of the resolution used 
or example the step size m the pressure intervals is 100 psi. However, linear interpolation is 
used to calculate values not found in the tables which reduces the error. Table 1 lists the data for 
a ypical equilibrium point EQ { corresponding to a given set of outflow requirements. A Matiab 
program [7] calculates the- required valve positions. Later in this report we will simulate the model 
w en the objective of a controller is to steer the system from EQ { to EQ f in a smooth manner, 
e next step is to obtain the linearized model at a chosen equilibrium point. Consider the 

mear mo e (7) valid m the neighborhood of the equilibrium EQi listed in Table 1. The matrices 
A and B are evaluated to be 


-17.52 

-3393.08 1 

P — 

' 25.11 

332.74 

-42.01 1 

. -0.1597 

-43.45 J ’ 

JD — 

' 0.651 

0.6612 

-0.4919 j 


For an output y — w e . the linear approximation is 

Welin = C vfiow x -f Df iow u ve = [0.2386 66.84] x + [l.2Z}u ve 

In the figures that follow, we present the behavior of the mixer obtained from the nonlinear 
mo e . n the absence of real data, we shall refer to the nonlinear SIMULINK simulation as the 
real system. Several simulations were performed to evaluate the linear model. However, one has to 
eep m mind that the linear model is only valid in the vicinity of the point where the linearization 
was performed. In other words, the linear model’s approximation gets worse as the state deviates 



from the linearization point, and it in actnally an invalid representation of the true behavior. The 
linear model will also be used later in the report to design linear controllers. 

• Disturbance . Effect - No Control. We ran the model under a 20 percent disturbance in the 
gas and liquid valves that take the state away from the initial equilibrium point. When the 
disturbance disappears, the mixer returns to the equilibrium. It was observed that the linear 
model captures the important dynamic characteristics of the system very accurately. 

• Tracking and Disturbance Rejection - Open Loop Control. The behavior of the mixer model 
was evaluated as the valves are steered in an open-loop fashion from the initial to the final 
positions. These were generated to simulate the manner in which an operator might open/close 
the valves m response to a visual monitoring of pressure, temperature, or flow sensors. The 
simulations included a 20 percent perturbations on the liquid and gas valve positions. This 

simulation run gave us a clear reminder that the linear model is only good in the vicinity of 
the linearization point. 


4 Point-to-Point Tracking Control Design 

As discussed in [6j, the mixer-valve system is a multi-input, multi-output system thus having a 
transfer function for each input-output pair. It is in general difficult to design PID-style controllers 

for each loop because of the interactions between channels. The approach chosen in this work is 
based on the stat 


l_ C i . n r+a m _ 


atner than on the individual transfer functions. 


4.1 Single Output Tracking 

Consider the problem of exit flow tracking and define a tracking error 

e(£) = w er - w e (t) ^ 

where the reference or desired constant exit flow is w„. We seek to design a controller that steers 

the error g.gnal signal to zero, that is, eft) - 0 as i - os. To that end, define an augmented state . 
vector 


X = 


x 

e 


that satisfies the differential equation 

a o ] r 

^ n x + 

'-'flow 0 0 0 — D 


X = 


flow J 


it. = AX + Bit 


( 10 ) 


( 11 ) 


Given the system (11), a performance criterion of the form 

roo , 

J ~ J Q dt 



is minimized. For technical reasons, the matrix Q must be positive-semidefinite and penalizes how 
far the state X[t) deviates from the origin. Similarly, the matrix R must be positive-definite and 

penalizes how much control efibri is used. The resulting optimal controller is given by the linear 
full-state feedback law 

■ - - KX[t) (12) 

where the matrix K is constant and is found a-priori by solving the steady-state Riccati Equation 
A T P + PA- PBR~ 1 B t P + Q = 0 =£■ K = R~ l B T P 

The closed-loop system 
X = [A - BK) X 


is guaranteed to be asymptotically stable, that is, X(t) 0, t -* co. Note that driving X to zero 
is equivalent to driving the error e to zero as required, and the mixer state x to a constant. Finally, 
integrating both sides of (12) leads to the optimal multi-input/single-output proportional-integral 
(MISO-PI) controller expression 


u(t) = — Kpx(t) 


K, 


Jo 


\w P 


w, 


■Ml dr 


(13) 


For example, selecting the following weighting matrices: 


1 

0 

0 


' 0.1 

0 

0 

0 

1 

0 

; R = 

0 

1 

0 

0 

0 

20 


0 

0 

0.1 


so that it is equally important to steer x x and x 2 to zero, but it is 20 times more important to drive 
the exit flow error to zero. Similarly, the liquid and exit valves have 10 times more authority than 
the gas valve in achieving the stated goal. The optimal gain matrix is found to be 



0.4679 

22,5556 

-12.1274 ' 


K - 

0.8992 

-9.3271 

-0.2184 | 


. 

-0.8653 

13.2205 

-7.2421 J 


which 






0.4679 

22.5556 ' 


’ -12.1274 ' 

K P = 

0.8992 

-9.3271 

; Kj 

= 

-0.2184 


_ -0.8653 

13.2205 


-7.2421 


The final equilibrium point was selected to be that corresponding to a 20 percent increase in exit 
flow, from w e = 37 Ibm/s to w e — 44.4 Ibm/ s. The resulting valve, mixer states, flows, and linear 
model performance was evaluated by simulations. There are (small) step changes that appear in 
the temperature and pressure plots caused by the resolution of the thermodynamic table and not 
by the controller. The PI controller successfully takes the exit flow to its final required value. 



Simulation results under the same PI controller and with the added perturbations on the liquid 
and gas valves show t hat the c ontroller recovers from these perturb ations . This b eha vior co uld not 
be achieved under open-loop schemes. However, the mixer states, and the valves do not go to their 
unique, pre-calcu-lated final values as would be given by a steady-state equilibrium calculation. This 
is because the controller is designed to track the exit flow step reference only. This is corrected in 
the next section. 


4.2 Multi Output Tracking 


Next, we design an optimal multi-input/multi-output PI (MIMO-PI) controller that simultaneously 
regulates the exit flow, the mixer pressure, and the exit temperature to arbitrary constant set points 
? Pvt ■} T er , respective!}-. The new augmented vector is 



x 

w er ~w e (t) 

Pvt ~ P v {t) 

T er - T e (t ) 


that satisfies the differential equation 



- 4 

0 

0 

0 



B 

■ 

X = 

C flow 

0 

0 

0 

X -r 

0 

0 

Dflow 


-Cpa 

0 

0 

0 


0 

0 

0 


~Cte 

0 

0 



0 

0 

0 . 


Minimizing the quadratic cost introduced in the previous section leads to the linear MI MO - PI 
controller 


u{t) = 


— Kpxit ) 



w eT - w e (r ) 
Pvt - P v (r) 
Ter - T e (r) 


dr 


As a simulation example, select the following weighting matrices; 


Q - 


1 0 0 0 0 
0 10 0 0 
0 0 20 0 0 
0 0 0 1 0 
.0 0 0 0 200 


R = 


1 0 0 
0 1 0 
0 0 1 


( 15 ) 


so that it is equally important to steer aq, aq, and the mixer pressure error to zero, it is 20 times 
more important to drive the exit flow error to zero, and it is 200 times more important to drive 
the exit temperature error to zero. Several simulations indicated that the temperature tracking is 



tie slowest of the three and hence the more stringent requirement. This time, all valves have equal 
authority m achieving the stated goals. The resulting optimal gain matrk is found to be 



-0.0719 

0.9454 

-0.0439 


46.93 

-9.8472 

-11.455 


-1.8728 

-0.384 

-4.0429 


-0.8897 

-0.1607 

0.4274 


2.5733 

-13.91 

0.1296 


The set point requirements are a 20 percent increase in exit flow from 37 to 44.4 Ibm/s ; a 100 psi 
increase m mixer pressure from 6000 psi ; and a 2 °F increase in exit temperature from -240 °F. 
These are introduced as step commands at time t = 0.3 sec. The valve disturbances are set on 
and a rate limiter with a slew rate at 50 is put in series with each valve to include another level 
of realism m the simulations. The simulation plots are in Figure 2 through 5. The effect of the 
rate limiter is an increase in overshoots and reduced overall performance quality while precluding 
unrealistically fast valve commands rises especially at the start of the command. The outputs are 
controlled as required. It is expected that real valves will not only exhibit rate limiting but other 
linear and nonlinear characteristics. These will be described in future reports. 


5 Conclusions and Further Research 

An analytical small-signal model of the high-pressure H 2 mixer and control valves has been derived. 
Tne model is suitable for a variety of analysis and control designs valid in the vicinity of an operating 
or equilibrium point. A classical approach to control design would involve the consideration of a total 
of nine single-input, single-output transfer functions. Unless the control loops are weakly coupled, 
spch a design is difficult to complete because of the interaction among the control ch ann els Instead, 
a state-variable control design has been pursued that results in an optimal multi-input/multi-output 
proportional/mtegral law that can be implemented in real-time. As illustrated in the report, it is 
straightforward to design the controller and it is possible to assign different authorities to each 
control valve. Precise and simultaneous tracking of constant setpoint changes in exit temperature, 

exit flow, and mixer pressure were achieved in the presence of constant and unmeasurable additive . 
valve position perturbations. 

Topics for further research include: 

• Introduce, a valve model that exhibits valve dynamics, and other nonlinearities such as friction 

and backlash. The linear model has to be recomputed because the valve introduces additional 
states. 

. Investigate the combined tank/mixer control problem. This is a six-state and four control-valve 
problem. Clearly, the quality of the out-flow properties is affected by the precise regulation of 
both the tank and mixer internal dynamic states. 

• Use heuristics-based set-point relegation to determine control valve authority. 

• Last but definitely not least, we must validate these findings with real data. 
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iaDle 1: Initial Equilibrium Data. Underlined Values are Computed. (IN=Into Valve: OUT=Out of 
Valve), zi in BTU/lbm: z 2 in lbm/ft 3 . 
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Appendix A 
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To determine the linear perturbation model (7), it is necessary to find the analytic expressions 
o t e indicated partial derivatives. For simplicity we include only those terms that do not evaluate 
to zero at the equilibrium point. Omitting the details, these are found to be: 
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Figure 2: MIMO PI Control vs Time (sec) (with valve 
State Trajectory: Start at o; End at x. 




perturbations). Bottom right plot is the resulting 





Figure 3: Mixer Behavior vs Time (sec). MIMO PI Control (with valve perturbations) 








